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Abstract: 

Heavy quark parton distribution functions (PDFs) play an important role in several 
Standard Model and New Physics processes. Most analyses rely on the assumption that the 
charm and bottom PDFs are generated perturbatively by gluon splitting and do not involve 
any non-perturbative degrees of freedom. It is clearly necessary to test this hypothesis with 
suitable QCD processes. Conversely, a non-perturbative, intrinsic heavy quark parton 
distribution has been predicted in the literature. We demonstrate that to a very good 
approximation the scale-evolution of the intrinsic heavy quark content of the nucleon is 
governed by non-singlet evolution equations. This allows us to analyze the intrinsic heavy 
quark distributions without having to resort to a full-fledged global analysis of parton 
distribution functions. We exploit this freedom to model intrinsic bottom distributions 
which are so far missing in the literature in order to estimate the impact of this non- 
perturbative contribution to the bottom-quark PDF, and on parton-parton luminosities 
at the LHC. This technique can be applied to the case of intrinsic charm, albeit within the 
limitations outlined in the following. 

Keywords: Parton distribution functions, PDFs, heavy quark PDFs, intrinsic charm, 
intrinsic bottom, LHC 
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1 Introduction 

Heavy quark parton distribution functions (PDFs) play an important role in several Stan¬ 
dard Model (SM) and New Physics (NP) processes at the CERN Large Hadron Collider 
(LHC). In particular, several key processes involve the bottom quark PDF, e.g. tW, tH^ 
production, associated b plus WjZjH boson production or Hbb production [1]. In the 
standard approach employed by almost all global analyses of PDFs, the heavy quark distri¬ 
butions are generated radiatively, according to DGLAP evolution equations [2-4], starting 
with a perturbatively calculable boundary condition [5, 6] at a scale of the order of the 
heavy quark mass.^ In other words, there are no free fit parameters associated with the 
heavy quark distribution and it is entirely related to the gluon PDF at the scale of the 
boundary condition. As a consequence, the uncertainties for the heavy quark and gluon 
distributions are strongly correlated; this has been discussed in the context of inclusive 
Higgs production at the Tevatron and LHC [12]. 

However, a purely perturbative, extrinsic, treatment where the heavy quarks are ra¬ 
diatively generated might not be adequate; in particular for the charm quark with a mass 
rric — 1.3 GeV which is not much bigger than typical hadronic scales but also for the bot¬ 
tom quark with a mass rrn, ~ 4.5 GeV. Indeed, there are a number of models that postulate 

^The most common approach is to use one of the general mass variable flavor number schemes (GM 
VFNS) such as ACOT [7, 8] or FONLL [9, 10]. For a review of the treatment of heavy quarks in PDF 
global analyses see e.g. [11]. 
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a non-perturbative, intrinsic, heavy quark component which is present even for scales Q 
below the heavy quark mass m. In particular, light-cone models predict a non-perturbative 
(’intrinsic’) heavy quark component in the proton wave-function [13, 14] and similar ex¬ 
pectations result from meson cloud models [15-17]; an overview of different models can be 
found, e.g., in [18]. Predictions of these models together with the EMC charm data [19, 20] 
motivated first theoretical analyses of the intrinsic charm (IC) content of the proton [21- 
23]. These first analyses were not global as they concentrated only on the possibility of 
explaining the EMC data. Later the CTEQ collaboration performed the first fully global 
analyses of PDFs including the IC possibility [24, 25]. These studies gave the first esti¬ 
mate, based on an array of different data sets, of how big the intrinsic charm could be. 
The possibility of IC was also considered by the MSTW group [26] . Most recently, two new 
global PDF analyses dedicated to IC have been performed:^ (i) the CTEQ collaboration 
has updated their previous work [28] using the CTIO NNLO framework [29], and (ii) an 
analysis of Jimenez-Delgado et al. [30] which is interesting as it uses less strict kinematic 
cuts which allow for the inclusion of low-Q, high-x data that should be more sensitive to 
the light-cone inspired IC component.^ These two most recent analyses set significantly 
different limits on the allowed IC contribution, partly because of the very different tolerance 
criteria which are used to define the range of acceptable fits. These differences highlight 
the utility of the techniques discussed in this paper as we can freely adjust the amount of 
IC/IB contributions (within limits, which we will quantify) without having to regenerate 
a complete global analysis for each case. 

It is essential to experimentally test the heavy quark PDFs, both the extrinsic and 
intrinsic components. One observable which is directly sensitive to an IC component is 
the deep inelastic charm structure function F^ixjQ'^). So far, the EMC data for F 2 is 
the only measurement of the charm structure function in the relevant (x, Q^) region which 
is sensitive to a large-x IC component; this is the only DIS data cited as evidence for 
intrinsic charm. The HERA data on and exist only for x < 0.1 and provide no 
constraints on a large-x IC or IB; only the inclusive structure function T 2 (x, Q^) at large- 
X measured quite precisely at HERA provides limited information on the large-x charm 
PDF. At a future Electron-Ion Collider (EIC) the large-x charm structure function would 
be accessible, and the rate for charm production at x > 0.1 could be increased by up to 
an order of magnitude due to the presence of a large-x non-perturbative intrinsic charm 
component in the nucleon [32] . Alternatively, the IC could be searched for by measuring the 
Callan-Gross ratio R{x,Q^) = Fl/Ft [33], or by studying angular distributions [34, 35]. 

In hadronic collisions, a promising way to constrain models of IC is the measurement of 
inclusive charm hadron production , D~^, , Ac, ■ ■.). Predictions for such processes 

were obtained in the general-mass variable-flavor-number scheme (CM-YENS) [36-38] for 

^From private communication we also know that the NNPDF collaboration will soon release an IC 
analysis based on the NNPDF3.0 framework [27]. 

® Of course this requires including additional corrections (target mass and higher twist) as the leading 
twist approximation does not necessarily hold in this kinematic regime. Also, as some of the data used is on 
heavy nuclear targets, nuclear corrections were employed as well. For an instructive discussion concerning 
this work see [31]. 
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the Tevatron at a center-of-mass (cms) energy of 1960 GeV, for the Relativistic Heavy 
Ion Collider (RHIC) at cms energies of 200 GeV (RHIC200) and 500 GeV (RHIC500), 
and for the LHC at a cms energy of 7 TeV (LHC7) [39-42], The IC charm effects can 
be particularly large at the RHIC200 and at the LHC7 at forward rapidities where the 
differential cross section can be enhanced by a factor of up to 5 compared to the prediction 
with a radiatively generated charm PDF [40] . Recent results from LHCb [43] can serve as 
an example of such data which are, however, not yet precise enough and do not extend 
to sufficiently large transverse momenta at the most forward rapidities (4 < y < 4.5) 
to be conclusive. Another process which is very sensitive to the heavy quark PDF is 
direct photon production in association with a heavy quark jet [44]. Data from the DO 
experiment at the Tevatron [45, 46] overshoot the standard NLO QGD predictions [47] at 
large transverse photon momenta; the inclusion of an intrinsic heavy quark component in 
the nucleon can reduce the difference between data and theory, but not fully resolve it. In 
fact, at the Tevatron the qq channel becomes important at large transverse momenta, and 
higher-order corrections to this channel not included in Ref. [47] might explain (part of) 
this discrepancy. The gg-channel does not play an important role at pp colliders; therefore, 
measurements of this process at RHIG and the LHC probe the heavy quark PDFs in 
different regions of the momentum fraction x and could shed more light on the current 
situation."^ A detailed study of 7 -|- Q production at the LHC operating at '/S = 8 TeV 
(LHC8) was performed in Refs. [50, 51]. There it was shown that the existence of IC in 
the proton can be visible at large transverse momenta of the photons and heavy quark jets 
at rapidities 1.5 < jy-yj < 2.4, jyd < 2.4. Indeed, for the BHPS model [13] the cross section 
can be enhanced by a factor of 2 or 3 for p7p > 200 GeV. However, the cross section is 
already quite small in this kinematic range so that this measurement will be statistically 
limited. The ideal place to observe or constrain intrinsic charm would be A Fixed Target 
ExpeRiment using the LHC beams (AFTER@LHC) [52-55] due to the lower cms energy 
{y/S = 115 GeV) together with a very high luminosity; for a review see [56]. 

While there are at least a few global analyses which allow for an intrinsic charm 
component in the nucleon [24, 25, 28, 30], studies of intrinsic bottom PDFs have not been 
performed at all. The main purpose of this paper is to outline a technique which can 
provide IB PDFs for any generic non-IB PDF set; we can then directly compare the IB 
PDFs with the non-IB PDFs to gauge the impact of the non-perturbative IB component of 
the nucleon structure on 6-quark initiated processes. Our approach exploits the fact that 
the intrinsic bottom PDF evolves (to an excellent precision) according to a standalone 
non-singlet evolution equation. Furthermore, due to the small momentum fraction carried 
by the IB PDF, the evolution of the other partons is essentially not disturbed by the IB 
component. These two observation allows us to compute the IB PDF without the need to 
perform a complete global analysis of PDFs. Thus, we can easily obtain a matched set of 
IB and non-IB PDFs. Note that because existing data entering global analyses of proton 
PDFs do not constrain the IB PDF, it would not be useful to try and obtain information 

^In addition, the measurements at the LHC would provide an important baseline for j + Q production 
in pA [48] and AA collisions [49] . 
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on the IB content of the nucleon using a global fit. 

The rest of this paper is organized as follows. In Sec. 2, we demonstrate that the 
scale-evolution of the intrinsic PDF is governed by a non-singlet evolution equation. We 
then propose suitable boundary conditions and perform a number of numerical tests of 
the quality of our approximations. In particular, we investigate to what degree the gluon 
distribution is perturbed by the presence of the intrinsic component. In Sec. 3, we use 
the IB PDFs to obtain predictions for parton-parton luminosities relevant at the LHC. In 
Sec. 4, we discuss our results and assess the impact of the intrinsic component. Finally, in 
Sec. 5, we summarize our results and present conclusions. 

2 Intrinsic heavy quark PDFs 
2.1 Definition 

In the context of a global analysis of PDFs the different parton flavors are specified via a 
boundary condition at the input scale //q which is typically of the order 0(1 GeV). Solving 
the DGLAP evolution equations with these boundary conditions allows us to determine the 
PDFs at higher scales > fig- The boundary conditions for the up, down, strange quarks 
and gluons are not perturbatively calculable and have to be determined from experimental 
data. From this perspective, it is meaningless to decompose the light quark and gluons 
PDFs into distinct (extrinsic and intrinsic) components. The situation is different for the 
heavy charm and bottom quarks; the boundary conditions have been calculated perturba¬ 
tively and resum to all orders collinear logarithms associated with the heavy quark lines at 
fixed-order in perturbation theory. A non-perturbative (intrinsic) heavy quark distribution 
Qi can then be defined at the input scale /tq as the difference of the full boundary con¬ 
dition for the heavy quark PDF Q and the perturbatively calculable (extrinsic) boundary 
condition Qq: 


Qi(x, Ho) := Q{x, Ho) - Qo(x, Ho), 


( 2 . 1 ) 


where Q = c or Q = b. At NLO in the MS scheme, the relation in Eq. (2.1) gets fur¬ 
ther simplified if the input scale ho is identified with the heavy quark mass mg because 
Qo(x,mQ) = 0 at NLO. In this case, any non-zero boundary condition Q(x,mQ) / 0 can 
be attributed to the intrinsic heavy quark component. This simplihcation, however, does 
not hold a priori for a different factorization scheme, nor at NNLO and beyond. 

2.2 Evolution 

In the following, we demonstrate that to a good approximation the intrinsic heavy quark 
distributions are governed by non-singlet evolution equations. Denoting the vector of 
light quarks as ‘g’ and the heavy quark distribution by ‘Q’ (where Q = c or Q = b) 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations read [2-4] 


Q = Pqq ®g + Pgq ®q + Pqq ® Q 


g = Pgg® g + Pgq® q + PgQ 0 Q 


q = Pqg ^ g + Pqq® q + PqQ ® Q 


( 2 . 2 ) 

(2.3) 

(2.4) 
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with the splitting functions PQg{x) = Pqg{x), Pqq{x) = Pqg{x), PQq{x) = Pq'q{x), etc. in 
the massless MS scheme which are known up to three-loop order [57, 58]. 

Next we substitute Q = Qo + Qi where Qo denotes the usual radiatively generated 
extrinsic heavy quark component and Qi is the non-perturbative intrinsic heavy quark 
distribution:^ 


g = Pgg® g P Pgq® q +PgQ® Qo +PgQ^SrQ{, (2.5) 

q = Pqg® g + Pqq® q + PqQ ® Qo P PgQ-®r^, (2.6) 

Qo P Qi = PQg ® g P PQq ® qP Pqq ® Qo P Pqq ® Ql ■ (2.7) 

Neglecting the crossed out terms which give a tiny contribution to the evolution of the 

gluon and light quark distributions the system of evolution equations can be separated 
into two independent parts. For the system of gluon, light quarks and extrinsic heavy 
quark {g,q,Qo) one recovers the same evolution equations as in the standard approach 
without an intrinsic heavy quark component 


g = Pgg®gPPgq®qPPgQ®Qo, (2.8) 

q = Pqg® g P Pqq® qP PqQ®Qo, (2.9) 

Qo = Pqg ® g P Pqq ® qP Pqq ® Qo ■ (2-10) 


For the intrinsic heavy quark distribution, Qi, one finds a standalone non-singlet evolution 
equation 

Qi = Pqq®Qi- (2-11) 

In a global analysis with intrinsic heavy quark PDFs, using the exact evolution equa¬ 
tions (2.2)-(2.4), the parton distributions satisfy the momentum sum rule 


dx x \ g P /^{qi P qi) P Qo P Qo P Ql P Ql ] — i • 


( 2 . 12 ) 


Allowing for a small violation of this sum rule it is possible to entirely decouple the analysis 
of the intrinsic heavy quark distribution from the rest of the system. The PDFs for the 
gluon, the light quarks and the extrinsic heavy quark can be taken from a global analysis in 
the standard approach using Eqs. (2.8) - (2.10) where they already saturate the momentum 
sum rule 


dx X I 5 -h y^iqi + gi) + Qo + Qo 1 = 1 ■ 


i 


(2.13) 


On top of these PDFs the intrinsic heavy quark PDF can be determined in a standalone 
analysis using the non-singlet evolution equation (2.11). This induces a violation of the 
momentum sum rule by the term 


dx X (Ql -b Ql) 


(2.14) 


^Strictly speaking, the decomposition of Q into Qo and Qi is defined at the input scale where the 
calculable boundary condition for Qo is known. Consequently, Qi Q — Qo is known as well. Only due 
to the approximations in Eqs. (2.5) and (2.6) it is possible to entirely decouple Qo from Qi so that the 
decomposition becomes meaningful at any scale. 
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which, however, is very small for bottom quarks.® We will perform numerical checks of the 
validity of our approximations in Sec. 2.5 after having discussed the boundary conditions 
for the intrinsic heavy quark distribution. 

2.3 Modeling the boundary condition 

The BHPS model [13] predicts the following x-dependence for the intrinsic charm (IC) 
parton distribution function; 

ci(x) = ci(x) oc x^[6x(l + x) Inx + (1 — x)(l + lOx + x^)] . (2-15) 

Conversely, the normalization and the precise energy scale of this distribution are not 
specified. In the CTEQ global analyses with intrinsic charm [24, 25] this functional form 
has been used as a boundary condition at the scale Q = rric leaving the normalization as a 
free fit parameter. 

We expect the x-shape of the intrinsic bottom distribution 6i(x) to be very similar 
to the one of the intrinsic charm distribution. Furthermore, the normalization of IB is 
expected to be parametrically suppressed with respect to IC by a factor m^/ml ~ 0.1. 
Therefore, because the scale of the boundary condition is not fixed, the following two 
ansatzes for bi can be considered 

9 

7TL 

Different Scales : 6i(x,mfc) = —|ci(x,mc), (2-16) 

ml 


Same Scales : 6i(x,mc) = — %c\{x,mc) ■ (2-17) 

ml 

In the following we use the Same Scales boundary condition, Eq. (2.17), which remains 
valid at any scale Q. 

In this case, since ci = c —cq, it is possible to construct the IB PDF from the difference 
of the CTEQ6.6c and the standard CTEQ6.6 charm PDFs at any scale without having to 
solve the non-singlet evolution equation for the IB PDF. We will compare the two boundary 
conditions in Eqs. (2.16) and (2.17) in Sec. 2.5. Finally, let us note that it would be no 
problem to work with asymmetric boundary conditions, ci(x) ci(x) and 6i(x) bi{x), 
as predicted for example by meson cloud models [16]. 

2.4 Intrinsic heavy quark PDFs from non-singlet evolution 

For the purpose of this analysis we used the approximation of Sec. 2.2 to produce standalone 
IC and IB PDFs that can be used together with any regular PDF set sharing the same 
values for the QCD parameters, such as the strong coupling or the quark masses. For the IC 
PDF we used Eq. (2.15) to define the initial x-dependence at the scale of the charm mass, 
and fixed the normalization to match the one predicted by the CTEQ6.6cO fit [25]. The IB 
PDF was generated using the Same Scales boundary conditions of Eq. (2.17) together with 
the same x-dependent input of Eq. (2.15). If not stated otherwise, the normalization for the 

®It is also acceptable in case of charm provided that the allowed normalization of IC is not too big. 
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Figure 1. Scale-evolution of the IC (left) and IB (right) PDF at NLO according to the non-singlet 
evolution equation (2.11) using the Same Scales boundary condition in Eq. (2.17) with rric = 1.3 
and mb = 4.5 GeV. Shown are results for = 1.69,10,100,1000, and 10000 GeV^. 


IB PDF was chosen to be identical to the IC case scaled down by a factor m^/ml = 0.083. 
Both PDFs were then evolved according to the non-singlet evolution equation (2.11) and 
the corresponding grids were produced.^ We show these distributions for selected values 
of the factorization scale (Q) in Fig. 1. As in our approximation, the evolution of the 
intrinsic charm and bottom PDFs is completely decoupled from the light quarks, gluons, 
and perturbative heavy quark components; the normalization of our PDFs can be easily 
changed by means of simple rescaling. However, for convenience we also produced a set 
with normalization corresponding to the CTEQG.Gcl fit [25] allowing for a larger intrinsic 
component. 

2.5 Numerical validation 

In order to test the ideas presented in Secs. 2.2 and 2.3, we use the CTEQG.Gc series of 
intrinsic charm fits which have been obtained in the framework of the CTEQ6.G global 
analysis [25]. The CTEQG.Gc series comprises 4 sets of PDFs including an intrinsic charm 
component. Two of them, CTEQG.GcO and CTEQG.Gcl, employ the BHPS model with 1% 
and 3.5% IC probability, respectively.® This corresponds to the values of 0.01 and 0.035 
of the first moment of the charm PDE, f dxc(x), calculated at the input scale Qo = = 

1.3 GeV. Eor a review of these models see Ref. [24]. In the rest of this article, we will 
follow the naming convention of the CTEQG.Gc hts in which a given fit is characterized by 
the value in percentage of the first moment of the charm distribution at the input scale, 

^The evolution was performed in Mellin-moment space using the PEGASUS package [59] at NLO in 
the MS scheme. Additionally the strong coupling was chosen so that it corresponds to the one used in the 
CTEQ6.6 set [25]. 

®The other two sets, CTEQ6.6c2 and CTEQ6.6c3, study a ’sea-like’ intrinsic charm with low and high 
strength, respectively. We don’t consider them in this paper because they are theoretically less motivated. 
Note also that the picture of a non-singlet intrinsic heavy quark distribution does not naturally apply, since 
these distributions are substantial in the small-* region. Additional numerical tests would be needed in 
this case. 
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(a) (b) 

Figure 2. (a) CTEQ6.6cO charm distribution function c{x, Q^) (solid lines) and the sum co(a;, <5^) + 
Ci{x, (triangles) where Cq is the radiatively generated CTEQ6.6 charm distribution and Ci is the 
non-singlet evolved IC using the BHPS boundary condition (2.15) with the same normalization as 
used for the CTEQ6.6cO charm distribution. Results are shown for the input scale = Qq = ra^ 
and the scale = 10000 GeV^. Fig. (b) shows the ratio of the curves in (a). 


e.g. 1% for CTEQ6.6cO. For convenience, we list below the first and second moments 
(calculated at the input scale) for the sets referred to in the following. 



Iq dx c{x) 

dx X [c(a:) + c{x)] =< x >c+c 

CTEQG.G 

0 

0 

CTEQG.GcO 

0.01 

0.0057 

CTEQG.Gcl 

0.035 

0.0200 


As can be seen the actual momentum carried by the charm in the CTEQG.GcO and CTEQG.Gcl 
fits is equal to ~ O.G% and 2% respectively. 

In the following we compare our approximate IC PDFs supplemented with the central 
CTEQG.6 fit, which has a radiatively generated charm distribution, with the CTEQG.GcO 
and CTEQG.Gcl sets where IC has been obtained from global analysis without the approx¬ 
imations of Sec. 2.2. 

In Fig. 2(a) the CTEQG.GcO charm distribution function c(x, Q^) is shown (solid lines) 
for two scales, = 1.G9 and 10000 GeV^, in dependence of x. The doted lines have 
been obtained as the sum of cq{x,Q‘^) -b ci{x,Q'^) where cq is the radiatively generated 
charm distribution using the CTEQG.G PDF and ci is the non-singlet evolved IC using 
the boundary condition (2.15) with the same normalization as used for the CTEQG.GcO 
charm distribution. As can be seen in the ratio plot. Fig. 2(b), the difference between 
the sum cq + ci and the CTEQG.GcO charm distribution is tiny at low Q^, and smaller 
than 5% at the higher Q^. In other words, the IC distribution ci evolved according to 
the decoupled non-singlet evolution equation is in very good agreement with the difference 
c — Co representing the IC component in the full global analysis. 
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Figure 3. (a) Comparison of the CTEQ6.6cO (solid line) and the CTEQ6.6 (triangles) gluon 
distributions, (b) Ratio of the CTEQ6.6cO and the CTEQ6.6 gluon distributions. The results are 
shown as function of x for two scales, = 1.69 and = 10000 GeV^. 


The inclusion of the intrinsic charm distribution will alter the other parton distribu¬ 
tions, most notably the gluon PDF. In order to gauge this effect we compare in Fig. 3 the 
gluon distribution from the CTEQ6.6cO analysis with the one from the standard CTEQ6.6 
fit. Fig. 3(a) shows the x-dependence of the gluon distribution for two scales, = 1.69 
and = 10000 GeV^; Fig. 3(b) shows the ratio of these curves. For small x (x < 0.1) 
the gluon PDF is not affected by the presence of a BHPS-like intrinsic charm component 
which is concentrated at large x. At x ~ 0.7, the CTEQ6.6cO gluon is suppressed by 
about 20% with respect to CTEQ6.6, and this is relatively insensitive to the value of Q^. 
We note that at large-x, the gluon distribution is already quite small and the uncertainty 
of the gluon PDF is sizable (of order of 40 - 50% for the CTEQ6.6 set). The difference 
between the gluon distributions is slightly enhanced when evolving from the input scale 
= 1.69 GeV^ to the electroweak scale = 10000 GeV^, but it is still much smaller 
than the PDF uncertainty. We conclude that for most applications, adding a standalone 
intrinsic charm distribution to an existing standard global analysis of PDFs is internally 
consistent and leads to only a small error. Moreover, for the case of intrinsic bottom which 
is additionally suppressed, the accuracy of the approximation will be even better. 

Another source of uncertainty is the choice of boundary conditions for the intrinsic 
distributions. For the IB PDF we have presented two equally compelling choices: the 
Different Scales ansatz of Eq. (2.16) and the Same Scales one of Eq. (2.17). In Fig. 4(a) 
we compare the Different Scales (dashed line) and the Same Scales boundary conditions 
(solid line) at the scale Q = mi, = 4.5 GeV. As expected, the Same Scales boundary 
condition leads to a softer distribution due to the evolution from rUc to The ratio of 
these two distributions varies between 2 at small x and ^ at large x (c/.. Fig. 4(b)).® Note 

®We show here the ratio for x > 0.1; at lower x values the intrinsic distribution bi is negligible compared 
to the perturbative bo component, as can be seen in Fig. 7. 


- 9 - 
































2.0 



Figure 4. (a) Comparison of Different Scales, Eq. (2.16), and Same Scales, Eq. (2.17), boundary 
conditions for the IB PDE at the scale Q = ruf, = 4.5 GeV. Fig. (b) shows the ratio of the two 
curves in (a). 


that the uncertainty due to the normalization is larger; for example, the ratio of the IC 
distributions ci{x,miy) from CTEQ6.6cl and CTEQG.GcO is about 3.5. Furthermore, there 
is also a freedom in the parametric m^/m^ factor depending on which values are used for 
the heavy quark masses. Therefore, in the following we use the Same Scales boundary 
condition for our numerical studies and consider the normalization as a free parameter. 

For completeness we also provide similar validation using the parton-parton luminosi¬ 


ties 


dCi. 

dr 


= 


dx 


1 + dij 


X L 


fi{x,n)fj{T/x,n) + fj{x,fL)fi{T/x,fl) 


(2.18) 


where r = xiX 2 - We can consider the production of a heavy final state particle H at the 
LHC {pp — )• H) of mass mjj = s, where s = X1X2S = tS so that r = rnjj/S. This will allow 
us to estimate the effects of our approximation on the physical observables as a function of 
the mass scale y/r = mn / We discuss the relation of the parton-parton luminosities 
to the actual cross-section in more detail in Sec. 3.2. 

In Fig. 5(a) we show the ratio of cc luminosities for the IC with two choices of nor¬ 
malizations; we compare the results obtained with our approach to the CTEQG.GcO and 
CTEQG.Gcl sets. We see that our error on the luminosities (the difference between corre¬ 
sponding solid and dashed lines) is smaller than 10% across the full range of r = mjj/S 
values. Additionally, this error is much smaller than the difference between CTEQG.G and 
scenarios with IC (distance between blue band and red/green lines). Similarly, in Fig. 5(b) 
we show the ratio of luminosities for the eg combination. In this case the error of our 
method is larger; it is around 10% for the IC 1% normalization, and it increases substan¬ 
tially for the case of IC 3.5% normalization. Note that for the IC 3.5% normalization, 
the deviation of our approximation from CTEQG.Gcl is of the same order as the CTEQG.G 
PDF error band. 
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Thus, we conclude our approach provides a good approximation for cc luminosities 
for an IC with either 1% or 3.5% normalization. For eg luminosities, this approximation 
is quite reasonable for an IC with 1% normalization; however, for 3.5% normalization, it 
is only sufficient to obtain a rough estimate of the effects. On the other hand, if the IC 
component is this large it should not be difficult to observe. 

Since the IB case can be obtained by scaling IC with the m^/m^ suppression factor, 
our approximation will work perfectly well for both the hh and the hg luminosities due to 
the smaller normalization. 





(a) (b) 

Figure 5. Parton-parton luminosities at the LHC with \/S = 14 TeV for (a) {c, c} and (b) {c,g} 
calculated using: our IC and CTEQ6.6 PDFs (solid lines), CTEQ6.6c PDFs (dashed lines). For 
comparison, we show the PDF uncertainty (blue band) for the CTEQ6.6 PDFs. The upper plots 
show ratio of the corresponding luminosities compared to the result obtained with the central fit of 
CTEQ6.6. The lower plots present the ratio TZ of luminosities calculated using our IC plus CTEQ6.6 
PDFs compared to the result with CTEQ6.6cO (CTEQ6.6cl) PDFs. Note that the red curves 
represent cases with 1% IC normalization and the green ones cases with 3.5% IC normalization. 


3 Possible effects of IC/IB on LHC observables 

In this section, we investigate the effects of intrinsic heavy quarks on observables at the 
LHC. We study the effects of both IC and IB on parton-parton luminosities at 14 TeV 
LHC. This allows us to assess the relevance of a non-perturbative heavy quark component 
for the production of new heavy particles coupling to the SM fermions. 

3.1 IB versus IC and prospects of observing IB 

The IB PDF is parametrically suppressed by a factor ml/mf ~ 1/10 compared to the 
IC PDF. However, whenever particles have couplings to the SM fermions proportional to 
the fermion mass, the m^/m^ suppression will be compensated by a factor of in the 
coupling. In addition, the radiatively generated (extrinsic) bottom PDF is typically a factor 
of ~l/5 to ~l/2 smaller than the extrinsic charm PDF as shown in Fig. 6. Combining 
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Figure 6. Ratio of the dynamically generated bottom and charm parton distributions at = 100 
and = 106 GeV^ as a function of x. The results for x > 0.7 are not reliable due to instabilities 
in the CTEQ PDF grids for the charm distribution (note however, that for that large x values the 
perturbative PDFs are nearly zero and are negligible). 




(a) (b) 


Figure 7. Ratio of intrinsic and dynamically generated heavy quark PDFs for (a) charm and (b) 
bottom quarks. The distributions Ci{x, Q"^) and 61 (a;, have been generated using the non-singlet 
evolution equation (2.11) with the BHPS form (2.15) (properly normalized to 1%) and the Same 
Scales boundary condition (2.17). co(c, Q^) is the CTEQ 6.6 charm distribution and bo{x, Q^) is the 
CTFQ 6 . 6 cO bottom distribution. As in Fig. 6 , the results in (a) are not reliable for x > 0.7 due to 
instabilities in the CTEQ 6.6 grids for very low values of the Cq distribution. 


these pieces together we find 61/60 < ci/cq which means that the possible effects due to 
the IB will be less pronounced than the ones due to the IC. The ratios ci/cq and 61/60 are 
shown for several scales in Figs. 7(a) and 7(b). 

These figures are useful to illustrate the impact of the intrinsic heavy quark distribu¬ 
tions on the physical observables. For this purpose, we define the ratio Kc = 1 + ci/cq 
and similarly Kb = 1 -|- 61/60 which measures the relative deviation expected due to the 
IC/IB components. For example, if the 6 -quark initiated subprocesses of an observable 
contributes a fraction rb to a cross section (say rb = 80%), the observable will be enhanced 
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by a factor [ri, x k;, + (1 — V},)] where Kb is evaluated at the x-value relevant for the specific 
process. 

We note that Kc is still sizable at Q = 100 GeV with Kc — 3 {kc — 8) for the case of 
1% (3.5%) normalization.^^ In contrast, the IB content is much smaller at this scale with 
Kb ^ 1.6 {Kb — 3.1) for the case of 1% (3.5%) normalization of the IC due to the 
factor. Therefore, we expect processes at the electroweak scale (or heavier scales) to be 
much less affected by the presence of an intrinsic bottom component. Nevertheless, in cases 
where a process is dominated by the 6-quark initiated subprocesses and where the large x 
region is probed (for example at large rapidities) an enhancement by a factor 1.6 (or even 
3.1) might be visible. In general, processes probing lower scales {Q < 20 GeV) and large 
X {x > 0.5) will be better suited to find/constrain the intrinsic bottom component of the 
nucleon. 


3.2 Parton—parton luminosities at the LHC 

We now turn our attention to the parton-parton luminosities and study the impact of a 
non-perturbative heavy quark component on these quantities for the LHG at 14 TeV. Using 
the factorization theorem of QCD for hadronic cross sections, one can express the inclusive 
cross section for the production of a heavy particle H as follows: 


(^pp^H+X 


E 


1 rl 


t/x\ 


dxidx 2 fi{xi, fi)fj{x 2 , 


(3.1) 


where r = xiX2 = mjj/S, S is the hadronic center of mass energy, and s = X1X2S is its 
partonic counterpart. fi{x,fi) denotes the PDF of parton i carrying momentum fraction x 
inside the proton. Finally, p is the factorization scale which in the following is identified 
with the partonic center of mass energy s = rnjp. Equation (3.1) can be re-written in the 
form of a convolution of partonic cross-sections and parton-parton luminosities [60, 61], 

\ dr ^ aij{s), (3.2) 


CTpp^H+X = / , 

b 


where ^) has been introduced in Eq. (2.18). All the results of this section have been 

obtained using the CTEQ6.6 PDF set [25] supplemented with our approximate IC and IB 
PDFs constructed using the procedure presented in Sec. 2. 

In Fig. 8(a) we show different parton-parton luminosities, dCij/dr, for the LHC at 
14 TeV (LHC14) as a function of -y/r = mn/VS- We choose the range of y/r to be 
[0.02,0.5] that corresponds to the production of a heavy particle of mass niH G [0.280, 7] 
TeV which is roughly the range of values that will likely be probed at the LHC14. As 
can be seen, at large ^/T, the parton-parton luminosities respect the following ordering: 
ug ^ uu > gg ^ go gb ^ cc > 66. Consequently, one can generally conclude that heavy 
quark initiated subprocesses play a minor role in most processes where a heavy state is 
produced. 


^°Note that we show only the case of 1% normalization, however, due to the scale-invariance, the 3.5% 
normalization can be obtained by applying a multiplicative factor 3.5. 
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Figure 8. (a) Different parton-parton luminosities as a function of ^fr = ranl^fS for the LHC 14 
calculated using CTEQ6.6 PDFs. For large r, the ordering of the curves is: ug ^ uu > gg ^ 
gc> gb'^ cc> bb. (b) Rescaled parton-parton luminosities {rnfdCij/dT) for the LHC14 calculated 
using CTEQ6.6 PDFs. For comparison, we also show the un-rescaled gluon-gluon luminosity. For 
large r, the ordering of the curves is: gg ^ gb > gc ^ bb > cc ^ ug ^ uu. Note that by coincidence 
the gluon-gluon luminosity, Lgg, agrees at the 10% level with the scaled gb luminosity, mlCgb, so 
that the two curves lie on top of each other in (b). 


One exception would be SM extensions where the couplings to the first two generations 
are suppressed or vanish so that the gb or bb channels can dominate; typically this is done 
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in order to avoid experimental constraints from low energy precision observables or flavor 
changing neutral currents. Of course, unless the couplings to the gb or bb channels are 
enhanced, these scenarios have tiny cross sections and will be difficult to measure at the 
LHC. 

However, if the couplings are enhanced by factors of the quark mass, the hierarchy 
of the contributions can change dramatically. This can happen when the heavy state 
has couplings to the Standard Model particles proportional to their masses such as the 
SM Higgs or the Higgs particles in 2HDM models. For example, in Fig. 8(a) we show 
the parton-parton luminosities with no enhancement factors; in Fig. 8(b) we show the 
same but with additional factors proportional to the heavy quark mass; the change is 
dramatic. Taking the quark masses into account, the high r region now exhibits the 
following hierarchy: gg gb > gc ^ bb > cc ^ ug ^ uu. In this case the heavy 
quark initiated subprocesses could play the dominant role, apart from the gg initiated 
subprocesses which would contribute via an effective, model-dependent, heavy quark loop- 
induced ggH coupling. 

To explore how the presence of IC and IB would affect physics observables with a 
non-negligible heavy quark initiated subprocesses, in Fig. 9, we present the parton-parton 
luminosities dC/dr, dC/dr with and without the intrinsic components. While the impact 
of the IC component for ^/r > 0.1 is clearly visible, the corresponding effect for the bottom- 
quark is smaller and lies inside the PDF uncertainty band. Generally, the enhancement is 
larger for the processes initiated by two heavy quarks {cc, 66} since the intrinsic component 
is then “squared” in the luminosities. 



Figure 9. Parton-parton luminosities dCcg/dr, dLcc/dr (left) and dCbg/dr, dCbi/dr (right) at 
the LHC14 as a function of yT = mH/V~S. Shown are results without an intrinsic heavy quark 
component using the CTEQ6.6 PDFs (solid lines) including the band due to the PDF uncertainties 
and the same quantities calculated with addition of IC normalized to 1% and the corresponding IB 
component (dashed lines). 

In order to precisely quantify the impact of the intrinsic components in Figs. 10 and 11 
we show the ratios of luminosities for charm and bottom with and without an intrinsic 
contribution for 1% and 3.5% normalizations. Furthermore, since there are no experimental 
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Figure 10. Ratio of cc luminosities (left) and eg luminosities (right) at the LHC14 for charm- 
quark PDF sets with and without an intrinsic component as a function of i/r = mn/'/S. The 
ratio for the cc luminosity (solid, green line) in the left figure reaches values of 50 at y/r = 0.5. In 
addition to the curves with 1% normalization (red, dashed lines) we include the results for the 3.5% 
normalization (green, solid lines) which was found to be still compatible with the current data [25] . 


constraints on the IB normalization, in Fig. 11 we also include an extreme scenario where 
we remove the usual m^/m^ factor; thus, the first moment of the IB is 1% at the initial 
scale rric- 

For the 1% normalization the cc luminosity ratio grows as large as 7 or 8 for ^/T = 0.5, 
and for a 3.5% normalization it becomes extremely large and reaches values of up to 50. 
From these figures we can clearly see that the effect of the 3.5% IC is substantial and can 
affect observables sensitive to cc and eg channels. As expected, in the case of IB the effect 
is smaller but for the bb luminosity the IB with 3.5% normalization leads to a curve which 
lies clearly above the error band of the purely perturbative result. In the extreme scenario 
(which is not likely but by no means excluded) the IB component has a big effect on both 
the bb and bg channels. 

4 Discussion 

We have demonstrated that the scale evolution of intrinsic heavy quark distributions (both 
charm and bottom) is governed by a non-singlet evolution equation to a very good approx¬ 
imation. Furthermore, the small intrinsic heavy quark distribution does not significantly 
influence the other parton distributions or the sum rules of a global analysis. This observa¬ 
tion holds to a very good precision for the intrinsic bottom case, but also works reasonably 
well for the intrinsic charm case (if the momentum fraction is not too large). Therefore, it 
is possible to perform a standalone analysis of the intrinsic heavy quark distribution and 
to combine it with the PDFs of a standard global analysis with dynamically generated 
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Figure 11. Ratio of luminosities at the LHC14 for bottom-quark PDF sets with different normal¬ 
izations of the intrinsic bottom component. The plot has been truncated, and the bb luminosity in 
the extreme scenario reaches about 17 at ^ = 0.5. 

heavy quark distributions^^. Note, this allows us to use any general PDF set and generate 
a matched IC or IB component without a global ht re-analysis. 

Based on this observation we have modeled an intrinsic bottom distribution and dis¬ 
cussed its effect on the relevant parton-parton luminosities at the LHC14. As a general 
rule, the effects of IB are less pronounced than the ones from IC due to the expected 
m^/m^ ~ 0.1 suppression factor. For example, we see from Fig. 7 (1% normalization) that 
whereas = I + ci/cq ~ 3 at scales Q ~ 100 GeV and x G [0.4,0.8] the corresponding 
factor for IB is relatively small: k;, = 1 + bi/hQ < 1.4. For a 3.5% normalization, the curves 
in Fig. 7 would be scaled by a factor 3.5, such that Kb < 2.4. 

We then turned to a discussion of the parton-parton luminosities where the main 
results can be seen in Figs. 10 and 11. 

The cc luminosities are strongly enhanced for the IC with both 1% and 3.5% normal¬ 
ization as the heavy quark factors enter quadratically; but the effect is smaller for the eg 
luminosity as there is only one factor of the heavy quark PDF. However, the effect is still 
signiheant; the 1% IC lies at the edge of the PDF uncertainty band, and the 3.5% IC yields 
a factor of 

As expected, the enhancement of the bb luminosity is much smaller compared to the 
cc case. For the 1% IB the curve lies well within the PDF uncertainty band. For the 
3.5% IB the results lies above this uncertainty band and predicts an enhancement of a 
factor 4 at ^/T = 0.5. For the bg luminosity, both the 1% IB and the 3.5% IB curves lie 
within the PDF uncertainty band. This band is largely driven by the uncertainty of the 
gluon distribution and might shrink in the future such that the enhancement due to the 
IB could become significant. For illustration, we have also included results for the extreme 

^^Needless to say, that the intrinsic heavy quark distribution could also be used together with the PDFs 
in a fixed-fiavor-number scheme where no dynamically generated heavy quark distribution is present. 
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assumption that the probability (first moment) for IB is 1%. In this case the enhancement 
is sizable and reaches a factor 17 for the bb case and a factor 4 for the bg case at ^/T = 0.5. 

In view of this, we can address the impact of IB on heavy new physics and certain 
electroweak processes where the bg or the bb channel plays an important role. For the 1% 
and 3.5% IB the enhancement for the ft^f-initiated subprocesses would be hidden within the 
PDF uncertainties. The 66-initiated subprocesses could be significantly enhanced in the 
case of 3.5% IB. The effect would be more pronounced for heavier states where however, 
the bb luminosity is extremely small so that a measurement would be limited by statistics 
even for models with enhanced couplings to the b quark. All in all, we conclude that the 
IB will have limited impact on searches for heavy new physics at the LHC. 

5 Conclusions 

In this article, we presented a method to generate a matched IC/IB distributions for any 
PDF set without the need for a complete global re-analysis. This allows one to easily carry 
out a consistent analysis including intrinsic heavy quark effects. Because the evolution 
equation for the intrinsic heavy quarks decouples, we can freely adjust the normalization 
of the IC/IB PDFs. 

For the IB, our approximation holds to a very good precision. For the IC, the error 
increases (because the IC increases), yet our method is still useful. For an IC normaliza¬ 
tion of 1-2%, the error is less than the PDF uncertainties at the large-x where the IC is 
relevant. For a larger normalization, although the error may be the same order as the PDF 
uncertainties, the IC effects also grow and can be separately distinguished from the case 
without IC. In any case, the IC/IB represents a non-perturbative systematic effect which 
should be taken into account. 

The method presented here greatly simplifies our ability to search for, and place con¬ 
straints upon, intrinsic charm and bottom compoments of the nucleon. This technique will 
facilitate more precise predictions which may be observed at future facilities such as an 
Electron Ion Collider (EIC), the Large Hadron-Electron collider (LHeC), or AFTER@LHC. 

The PDF sets for intrinsic charm and intrinsic bottom discussed in this analysis (1% 
IC, 3.5% IC, 1% IB, 3.5% IB) are available from the authors upon request. 
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